Phonons from neutron powder diffraction 
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The spherically averaged structure function /5>(|q|) obtained from pulsed neutron powder diffrac- 
tion contains both elastic and inelastic scattering via an integral over energy. The Fourier transfor- 
mation of Sdql) to real space, as is done in the pair density function (PDF) analysis, regularizes 
the data, i.e. it accentuates the diffuse scattering. We present a technique which enables the extrac- 
tion of off-center (|q| 7^ 0) phonon information from powder diffraction experiments by comparing 
the experimental PDF with theoretical calculations based on standard interatomic potentials and 
the crystal symmetry. This procedure (dynamics from powder diffraction(DPD)) has been success- 
fully implemented for two systems, a simple metal, fee Ni, and an ionic crystal, CaF2. Although 
computationally intensive, this data analysis allows for a phonon based modeling of the PDF, and 
additionally provides off-center phonon information from powder neutron diffraction. 
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For a variety of physical questions it is significant to ob- 
tain information on off-center phonons in crystals. This is 
particularly important in complex materials where local 
effects modify the macroscopic properties, as for example 
in high temperature superconductors, colossal magneto- 
resistive materials, ferroelectrics, intermetallic alloys and 
many more. Until this study the only available method to 
obtain off-center phonon data has been inelastic neutron 
scattering, which is intensity limited, hence time consum- 
ing, and for detailed studies, relies on the availability of 
large single crystals (triple axis measurements) . Here we 
show how to obtain similar data from powder neutron 
diffraction. 

Historically, the purpose of powder (polycrystalline) 
neutron diffraction has been the exact determination 
of the average crystal structure using modern crystallo- 
graphic analysis techniques like Rietveld refinement, see 
e. g. [Q. Within this type of refinement only a limited 
range of the momentum transfer q is necessary. More 
recently, however, the availability of pulsed sources has 
made possible the measurement of iS(|q|) up to very large 
values of q. The PDF analysis 0] has made use of this 
additional information to investigate local atomic devia- 
tions from an average crystallographic structure. 

To model the peak positions in the PDF analysis 
one calculates interatomic distances from the crystallo- 
graphic unit cell. The peak shape is commonly fitted to 
the experimentally obtained PDF, p exp (r), using Gaus- 
sians. Since the experimental PDF contains additional 
information from the diffuse scattering, observed differ- 
ences have been successfully attributed to local structural 
deformations like e. g. polarons H. This type of model- 
ing of the PDF does not take into account the intrinsic 
peak widths caused by coherent excitations in solids like 
phonons or spin waves, and how the peak widths are 
modified by a finite momentum transfer cut-off. Other 
attempts to relate the measured PDF peak widths to the 
intrinsic phonon dynamics in a real space approach 0] 



were not intended to address the inverse problem, i. e. 
extracting the phonon parameters from the PDF. In a 
real space approach, it is considerably more difficult to 
take into account possible extinction rules for inelastic 
scattering caused by the point group symmetry of the 
reciprocal lattice, as has been recently shown Also 
in reciprocal space the instrumental resolution function 
is easier taken care of. In addition to information on lat- 
tice dynamics we provide an extended modeling of the 
PDF based on the estimated phonon parameters. 

Triple axis neutron scattering is the most direct way to 
measure phonons. Indirect measurements on powders are 
either based on a detailed analysis of the peak shape in 
iS(|q|), as is done in the analysis of thermal diffuse scat- 
tering (TDS) @, or via a reconstruction of S'Hql, lo) from 
time-of-flight (TOF) experiments, see e.g. 0. Inelastic 
experiments are intensity limited, and the TDS analysis 
relies on a very good representation of the background. 
In this letter, we show how to extract information on 
off-center phonons in a reliable way from S'(lql) deter- 
mined from TOF experiments. This is achieved via a 
parameterization of the phonon dispersion curves and 
eigenvectors using a suitable model (as is standard in the 
presentation of triple axis data, see e.g. MM), from which 
we can calculate a theoretical S'dqj). The data are then 
regularized by transforming to real space to obtain a 
theoretical PDF, p theo (r). A reverse Monte Carlo pro- 
cedure is then used to estimate the parameters from a 
comparison with the experimental data, p exp (r), and to 
give data driven error bars. 

The measured quantity is the powder (angular) aver- 
aged structure factor 5"(|q|) (normalized by the neutron 
scattering lengths). For simplicity we assume that the 
incoherent effects can be neglected, and decompose the 
total structure factor into elastic and one-phonon inelas- 
tic contributions: S^(q) « Sf a h stlc (q) + S^q). Note 
that ^"^(q) contains all dynamic information as an in- 
tegral over frequency fllQfl. The elastic part is given by 
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and the one-phonon inelastic part by 
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where i^-(q) = Ed ^fe(q- <T d( k )) e_H ' ti(q)+4q,d ' the sum 
over G runs over all vectors of the reciprocal lat- 
tice, k describes a Brillouin zone integral over iVk 
points, d runs over atoms in the cell and j counts the 
branches, bd is the average scattering length and Ma 
the mass of the d'th atom, vq the unit cell volume, and 
W(q) is the Debye-Waller factor defined as Wd(q) = 

TOEkeB^^W^cothlf^lk)) . Thephonons 
enter eqns (||) and (||) via the dispersion relations Wj (k) 
and the eigenvectors cr^(k). The powder averaged struc- 
ture function 5(|q|) = i ]j q | =const dO q S(q) must be ob- 
tained numerically, and is CPU time consuming for the 
inelastic part of eqn(^). To compare to the experimental 
situation it is necessary to convolute S(|q|) with a suit- 
ably chosen resolution function |j] which depends on the 
instrument. 

FIG. 1. The experimental p exp (r) compared with the the- 
oretical p °(r) calculated for the expectation values {/,;} re- 
sulting from the DPD procedure for Ni at room temperature. 
<r(r) are the experimental errors in the x 2 estimator (Ul). 
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The extraction of the lattice dynamics described by 
LUj (q) and er^ (q) using a comparison of TOF experimen- 
tal data with a theoretical *S'( | q| ) is difficult in reciprocal 
space because 5(|q|) is a superposition of elastic scatter- 
ing, i.e. the Bragg peaks, and inelastic scattering, both of 



which are vastly different in amplitude. A direct compar- 
ison in reciprocal space is further complicated by the un- 
avoidable instrument background which is superimposed 
on the data. Both problems can be greatly reduced by 
transforming Sdql) into the real space pair density func- 
tion p(r) via Fourier transform 



P{r) = Po + 
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dq[S(\q\)-l]qsm(qr). (3) 



Pulsed sources provide a high enough momentum trans- 
fer to reduce the truncation error in eqn (^) and in turn, 
allow for a proper normalization of the data. Due to the 
properties of a Fourier transform the peaks in p(r) are all 
of comparable height which eliminates the problem with 
the large amplitude fluctuations in S"(|q|), i.e. it regular- 
izes the data analysis problem. As a useful byproduct, all 
features in 5(|q|) with inverse length scales larger than 
— , i.e. most of the background, are transformed into 
distances in p(r), which are smaller than the shortest in- 
teratomic distance ao, and can therefore be neglected in 
the parameter estimation. 

FIG. 2. The phonon dispersion curves for room temper- 
ature Ni. The symbols are triple axis data from jDJ. The 
shaded bands are the one-sigma error intervals of the dis- 
persion curves obtained by our procedure. The green bands 
result from just reproducing the PDF, and the red bands are 
from constraining the reverse Monte Carlo to agree with the 
measured elastic constants within their error bars. 
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We use a ^-functional in the parameter estimation 
process defined by 
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where a(r) is the error in p exp (r), which can be assumed 
to be H. The variational parameters are the param- 
eters describing the lattice dynamics {/j}, which depend 
on a theoretical model (see below for concrete examples). 
r a , r;, define the range over which the PDFs are com- 
pared. We use a Monte Carlo procedure Q which gives 
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statistical estimators for the parameters fi, and also er- 
ror bars (A/i) 2 = ff — fi . The A/j are controlled by 
the quality of the phonon model (a bad model also giving 
a large % 2 , i.e. systematic error), and by the quality of 
the data, the statistical error. Due to the nature of our 
data, the experimental PDF, the dynamics is obtained by 
averaging over the whole reciprocal space, not just along 
symmetry directions, as is usually done, when triple axis 
data are fitted. 

FIG. 3. The same as Fig. 1 but for CaFi at room temper- 
ature. 
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The parameterization of the lattice dynamics involves 
the choice of a theoretical model. There exists a large va- 
riety of models which describe the measured dispersions 
from triple axis data §||. To validate our procedure 
we chose to perform the above analysis on two materials 
with very different models for the lattice dynamics. The 
first example is fee Ni, whose phonon dispersion can be 
described by a simple force constant model p3|Jl^ ]. In 
Fig. 1, the PDF determined from the diffraction data is 
compared with the PDF calculated from the estimated 
parameters. The agreement is excellent and indicates 
that this procedure does indeed reproduce both the peak 
position and peak shape of a PDF measurement. A com- 
parison of the force constants resulting from our analysis 
with existing triple axis data is shown in Table 1. Given 
the force constants and their standard deviations we can 
easily calculate the resulting phonon dispersion curves. 
In Fig. 2 we show those curves together with published 
triple axis data [jl3| . Again the agreement is excellent 
proving that we managed to extract the dynamical infor- 
mation contained in «S'( | q| ) . Adding more information, 
e.g elastic constants, to the functional (Q) reduces the 
errors in our parameter estimation. The additional con- 
straints lead to smaller error bars, and a slight improve- 
ment in the overall shape of the dispersions. 

As a second example we took CaFi, since it involves 
a fairly ionic system, which necessitates a different in- 
teratomic force model. Also the phonons are a lot more 
complicated due to the presence of optical branches. The 



phonons were parameterized by ionic shell models, as de- 
scribed in Elcombe & Pryor |ll|. This is a more chal- 
lenging problem as we now have a much higher dynamic 
range (phonon energies up to 14THz ), and the optical 
modes show an unusually large dispersion. The results 
of the DPD analysis are presented together with some 
available triple axis data in Table 2 for the force con- 
stants, and in Figs. 3 and 4 for the comparison of PDFs 
and dispersions, respectively. Again, the agreement is re- 
markable and shows that we can obtain reliable off-center 
phonon information even for non-trivial systems. 

FIG. 4. The phonon dispersion curves for CaF^ at room 
temperature. The symbols are measured data from [16]. The 
shaded bands are one-sigma intervals of the dispersion curves 
obtained by our procedure. . 
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In this letter we have shown how dynamical informa- 
tion present in 5(|q|) data can be extracted. The success 
of the present analysis relies on high quality data and 
on intelligent modeling of the lattice dynamics. Noisy 
data or bad normalization will result in large error bars. 
The models can easily be extended and modified if the 
error analysis indicates through an unusually large \ 2 ~ 
value that the model is incompatible to the data. Such 
extensions may include more shells, or a totally different 
modeling based on pseudopotential calculations. As ad- 
ditional constraints to the functional in eqn.(|j) we can 
use further data from zone center experiments, like elas- 
tic constants and Raman data. At very low temperatures 
where the one-phonon processes are not thermally acti- 
vated, the information is not in the data. If necessary, 
more sophisticated Placzek corrections and the incoher- 
ent scattering can be included in the analysis. Our results 
are very encouraging, and we will pursue this effort on 
more complicated materials, where triple axis data are 
sparse or non-existent. 
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TABLE I. Generalized force constants / ± A/ for Ni at 
room temperature (in 10 4 dyn/cm) from reverse Monte Carlo 
runs with (a) the elastic constants constraints, (b) without 
, and from a triple-axis measurement Ref. iQ. The errors 
in the force constants in Ref. Jl4| are estimated at about 
±0.04 x 10 4 dyn/cm. The notation follows Ref. Q. 



(a) 



(b) 



Ref. 1 14 1 



1XX 
1ZZ 
1XY 
2XX 
2YY 
3XX 
3YY 
3XZ 
3YZ 



1.755 ± 
-0.054 ± 
1.878 ± 
0.067 ± 
-0.011 ± 
0.071 ± 
0.031 ± 
0.045 ± 
-0.010 ± 



0.018 
0.019 
0.017 
0.023 
0.008 
0.007 
0.004 
0.001 
0.005 



1.683 ± 
-0.039 ± 
1.869 ± 
0.122 ± 
-0.085 ± 
0.102 ± 
0.044 ± 
0.027 ± 
-0.009 ± 



0.024 
0.019 
0.159 
0.046 
0.031 
0.013 
0.006 
0.020 
0.005 



1.7319 
-0.0436 
1.9100 
0.1044 
-0.0780 
0.0842 
0.0263 
0.0424 
-0.0109 



TABLE II. The parameters of the shell model used for 
CaF2. The notation follows Ref. |l^| and the ionic charge 
Z(Ca) was fixed at 2e (shell model III in Ref. pi). 



Units 



RMC of PDF 



Ref. [15] 



Si 

A 2 
B 2 
A 3 
B, 
Aa 
B 4 

a(Ca) 
d(Ca) 
a(F) 
d(F) 



7« 



16.00 
-1.84 
1.30 
0.09 
0.22 
-0.28 
-0.18 



e 

A 3 

e 

A 3 



±0.25 
±0.09 
±0.11 
±0.01 
±0.02 
±0.03 
±0.02 
0.062 ± 0.006 
1.66 ±0.16 
0.25 ±0.02 
0.45 ±0.04 
0.065 ± 0.006 



15.12 
-1.70 
1.30 
0.09 
0.23 
-0.28 
-0.18 
0.061 
1.63 
-0.25 
0.45 
0.063 



±0.17 

±0.08 

±0.11 

±0.03 

±0.17 

±0.06 

±0.04 

±0.014 

±0.11 

±0.03 

±0.04 

±0.015 
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